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ABSTRACT 

We present a new series of supernova neutrino light curves and spectra cal- 
culated by numerical simulations for a variety of progenitor stellar masses (13- 
5OM0) and metallicities (Z = 0.02 and 0.004), which would be useful for a 
broad range of supernova neutrino studies, e.g., simulations of future neutrino 
burst detection by underground detectors, or theoretical predictions for the relic 
supernova neutrino background. To follow the evolution from the onset of col- 
lapse to 20 s after the core bounce, we combine the results of neutrino-radiation 
hydrodynamic simulations for the early phase and quasi-static evolutionary cal- 
culations of neutrino diffusion for the late phase, with different values of shock 
revival time as a parameter that should depend on the still unknown explosion 
mechanism. We here describe the calculation methods and basic results including 
the dependence on progenitor models and the shock revival time. The neutrino 
data are publicly available electronically. 
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Introduction 



Supernova explosion is one of the most spectacular events in the universe. It is not only 
the death of massive stars but also an engine for the evolution of the galaxies. Apart from 
the classification according to the spectroscopy, it is thought that there are two mechanisms 
of the explosion. While type la supernova is caused by the thermonuclear explosion of 
white dwarfs, the other types of supernovae are driven by the core collapse of massive stars. 
Unfortunately, however, details of the supernova explosion remain unsolved. 

As for the collapse-driven supernova, the physics which makes the explosion is not well 
understood. While many numerical simulations for the gravi t ational collapse of massive stars 



have been done so far (e.g. I0ttll2009l : [Thielemann et al.ll201ll : iKotake et al.ll2012l : |Jankall2012 

for recent reviews), there is no consensus on the explosion mechanism of collapse-driven su- 
pernovae. However, a rough sketch of the scenario shown below is widely accepted. The core 
of massive star becomes gravitationally unstable at the end of evolution and starts to col- 
lapse. The collapse is bounced by the nuclear repulsion force and the shock wave is launched. 
A supernova explosion is observed when the shock wave successfully propagates up to the 
stellar surface and expels the envelope material. Finally the neutron star or black hole is 
formed as a remnant. However, the general difficulty in numerical simulations for successful 
explosions is that shock waves tend to stall before blowing out the stellar envelopes. A possi- 
ble mechanism to revive the stalled shock waves is energy input by strong neutrino radiation 
from the newborn neutron stars ("the delayed explosion scenario"), but it is still a matter 
of debate whether this single process is sufficient for successful explosions. An aspherical 
hydrodynamic turbulence, such as the convective instability and standing accretion shock 
instability, may help the shock revival. Other proposed mechanism or physical processes to 
account for explosions include QCD phase transition, acoustic wave and magnetic field. 

The collapse-driven supernova is important also as a target of neutrino astronomy. As 
is well known, supernova neutrinos from SN1987A h ave been d etected by the two water 
Cerenkov detectors, Kamiokande H (IHirata et al.lll987l ) and 1MB (IBionta et al.lll987l ). Since 
the Large Magellanic Cloud, where SN1987A appeared, is ~50 kpc away from the Earth and 
detectors at the time were not as large as the present-day experiments, the event number 
was small (eleven for Kamiokande II and eight for 1MB). Nevertheless, many theoretical 
studies tried to extr act various in formation about physics of supernovae and/or neutrinos 
from these data (e.g. lRaffeltll2012l for a recent review). If a supernova occurs now near the 
Galactic center (~10 kpc away from the Earth), about 10,000 events will be detected by 
SuperKamiokande, which is the largest, currently operating neutrino detector around MeV 
( iBurrows et al.lll992t iTotani et al.l 119981 ). Therefore details of a collapse-driven supernova 
such as the explosion mechanism, as well as the neutrino physics such as mass and its hier- 
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archy, mixing and oscillation, or any exotic physics, could be investigated by high-statistics 
light curve and spectrum of a future supernova neutrino burst event. 

Another potential opportunity to observe neutrinos from supernovae is the relic back- 
ground radiation, which gives valuable cl ues about e.g., the c osmic star formation history 
and / or stellar initial mass function (e.g. iTotani et al.l Il996l : lAndo fc Satol |2004J : iBeacom 
20101 ). To predict the flux and spectrum of the relic supernova neutrinos, one must model 
the integrated spectrum of supernova neutrinos for various ranges of progenitor stars. 

Considering these situations, it would obviously be useful to provide the results of 
supernova neutrino emission calculated by state-of-art numerical simulations, for various 
types of progenitor stars. Such a comprehensive database of theoretical predictions can be 
used as templates to simulate supernova neutrino detection of a future neutrino burst event, 
or to construct a realistic prediction for the relic neutrino background. Howeve r, most of 
numerical simulations so far followed only within <1 s (e.g. iMarek fc Jankall2009l ) after the 
core collapse because solving the neutrino transfer with hydrodynamics for a long term is 
a numerically tough problem. Such simulations are not satisfactory for the above purpose, 
because a significant fraction of the total gravitational energy is emitted as neutrinos after 
this early phase, with a typical decay time scale of ~10 s. Such a simulation result for a 



long term was shown in iTotani et al.l (Il998[ ) . Recently, iFischer et al.l (120101 12012|) showed 
long term neutrino signals for two models, an ordinary collapse- driven supernova and an 
electron-capture supernova, which is a subcategory of collapse-driven supernovae. On the 
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2OI2I ). However, there is no comprehensive data set yet for numerical long-term supernova 



neutrino signals from various types of progenitor stars. 

In this study, we construct a publicly available database for numerical simulations of 
supernova neutrino emission, including eight simulations for normal single stars in the initial 
mass range of 13-5OM0 for two different values of stellar metallicity {Z = 0.02 and 0.004, i.e., 
the solar abundance and its 1/5, respectively). To follow the long-term evolution for many 
simulations, we use two different simulation methods: general relativistic neutrino-radiation 
hydrodynamic (z/RHD) simulations for the early phase, and general relativistic quasi-static 
evolutionary calculations of neutrino diffusion in a nascent neutron star for the late phase. 
Though the i/RHD simulations do not lead to a natural supernova explosion, we introduce 
a parameter, shock revival time, which reflects the unknown explosion mechanism, and 
phenomenologically connect the early and late phases by this parameter based on physical 
considerations. Although the connection between the two phases is not perfectly consistent 
as a physical simulation, this approach allows us to provide reasonably realistic supernova 
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neutrino light curves and spectra from the onset of a core collapse to ~10 s after that, for 
various progenitor stellar models. 

The organization of this paper is as follows. In Section |51 the progenitor models used 
in this study are described. The early phase of supernovae, from the onset of collapse to 
the shock propagation, is described in Section |3l In Section HI the quasi-static evolutionary 
calculations of neutrino diffusion in the dense core for the late phase are described. The 
connection between the early and late phases are discussed in Section O Section E] is devoted 
to a summary and discussion. 



Massive stars with the initial mass Mi^^t > IOMq are known to undergo gravitational 
collapse at the end of their lives. They synthesize heavy elements through some nuclear 
burning stages during the quasi-static evolutions. Finally the iron core, where all nuclear 
reactions become equilibrium to inverse reactions of themselves and various elements called 
"iron" (Fe, Ni, Co, Mn and so on) are created, is formed at the center. The iron core has den- 
sity of ^10^° g/cm^ and is supported by the degenerate pressure of relativistic electrons. On 
the other hand, the condition for the gravitational instability is written, using the adiabatic 
index, 7, as 



where p, p and s are the pressure, density and entropy, respectively. Since the adiabatic 
index of relativistic ideal gas is 4/3, the iron core is marginally stable. 

There are two processes which destabilize the iron core. The first one is an electron 
capture by protons belonging to nuclei: 



This reaction is caused because the sum of rest mass and kinetic energy of the relativis- 
tic electrons exceeds the mass difference between neutron and proton due to high density. 
Diminution of electrons cause a deficit in pressure and makes the iron core unstable. The 
second is the photodisintegration reaction of nucleus such as 



2. Progenitors of collapse-driven supernova 






(2) 




(3) 



This process works for high temperature as T ~ 5 x 10^ K. While the energy per nucleon, 
e, is minimum for ^^Fe, the Helmholtz free energy per nucleon, /, is minimized in finite 
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temperature. Because of / = e — Ts, matter becomes thermodjTiamically stable by the 
increase in entropy breaking nuclei. Nevertheless, this is an endothermic reaction so that 
the pressure does not increase so much as to suspend the core collapse. This is an onset of 
the collapse-driven supernova. 

The lower mass limit for the progenitors of collapse-driven supernova is determined 
whether the iron core is formed or not. Note that, for somewhat low mass cases, the core col- 
lapse may be triggered only by electron captures before Ne ignition, which is called electron- 
capture supernovae. They can be regarded as a subcategory of co llapse- driven supernovae 
because their scenarios after the collapse are thought to be similar. iPoelarends et al.l (|2008[ ) 
investigated stellar evolution sequences with initial masses between 6.5 and 13Mq and the 
solar metallicity using three differe nt codes and f ound that the lower mass limit for core col- 
lapse is 9-12Mq. Observationally, ISmartt et al.l (120091 ) found that it converges to 8 ± IMq 
investigating 92 samples which complete the collapse driven supernovae occurred in a fixed 
10.5-year period within a distance of 28 Mpc. 

The upper mass limit is determined whether the supernova explosion succeeds or not. 
It is theoretically uncertain because the precise understanding on the explosion mechanism 
is still absent. Nevertheless, generally to say, the density profile of progenitor is important 
because the shock wave should run through the star. What makes a situation more awkward 
is uncertainty in the theory of stellar evolution, the mass loss rate and convection. The 



effects of rotation and binary interaction may also affect the stellar evolution. iFryerl (119991 ) 
showed numerically that nonrotating stars with initial masses >40Mq fail to explode and 
form black holes directly while the mass loss was not considered in the progenitor models 
adopted by them. Note that, the black holes could b e also formed after the explosion, 
which is probably weak, via fallback accretion. Recently, lO'Connor &: OttI (120111 ) performed 
the core collapse simulations for many progenitor sets with the spherically symmetric (ID) 
models involving simplified neutrino transfer. They implied that the upper mass limit for 
the explosion depends severely on the evolutionary calculations, while their estimation was 
based on the calculations with an artificially increased energy deposition. 



Observational constraints on the upper mass limit are also obscure. ISmartt et al.l (120091 ) 
suggested that the maximum mass for the progenitors of type 11-P supernovae is ~2OM0 and 
the majority of more massive stars may col l apse q uietly to black holes where the explosions 
remain undetected. However, ISmith et al.l (120111 ). who investigated the observed fractions 
of the collapse-driven supernovae, claimed that they will produce other types of supernovae 
which are not II- P. On the other hand, massive stars with the mass of >35M0 lose the outer 
layer due to their strong stella r winds and form Wolf-Rayet stars, where the helium cores are 
exposed (IWoosley et al.ll2002l ). They are thought as progenitors of type lb and Ic supernovae 
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and, according to the light curve models by iNomoto et al.l (|2006[ ). the progenitors of some 
luminous type Ic supernovae called hypernovae have the masses of ^40Mq. They implied 
that the fate of massive stars beyond 25M(T) depends on t he stellar rotation and the hypernova 
progenitors are rotating. Note that, ISmith et al.l ( 120111 ) suggested that a certain fraction of 
such massive stars produce type Iln supernovae from luminous blue variables. They also 
pointed out the importance of binary interactions to account for type lb and Ic supernovae. 
Recently, extremely luminous supernovae such as SN2006gy (type Iln) and SN2007bi (type 
Ic) are observed and they are suggested to have a progenitor mass of >1OOM0, whereas they 
may be pair-instability supernovae, whose explosion mechanism is different from that of the 
collapse- driven supernovae. It should be emphasized that, as described above, the fate of 
massive stars is a hot subject under active discussion, but we deal progenitors provided from 
a single treatment for the stellar evolution in below. 

In this study, we prepare eight progenitor models with the initial masses Minit = ISM©, 
2OM0, 3OM0 and 5OM0, and the metallicities of Z = 0.02 (solar) and 0.004 (Small Mag- 
ellanic Cloud). They are computed by a Henyey type stellar evolution code wh i ch is fully 
coupled to a nuclear reaction network. This code is also used in lUmeda fc Nomotd (120081) and 



its descriptions, such as the treatment of convection, are given in 



for the mass loss rate, we adopt the sar ae model with "Case A" of lYoshida fc Umedal (1201 ll ) 



for the main sequence st age (see also. Ide Jager et al.l llQSSi : IVink et al 



empirical mass loss ra te (|de Jager et a 



the Wolf-Rayet stage (IKudritzki et al 



Umeda et al. 



20 



2h. As 



20011). We set the 
19881) scaled with the metallicity as (Z/0.02)°-^ for 



1989h . The effects of rotation and binary interaction 



are not taken into account. 

The density profiles of our progenitor models are shown in Figure [1] and the composition 
profiles for two reference cases, (Minit, ^) = (13M0,O.O2) and (SOM©, 0.004), are shown in 
Figure |21 Note that, as shown in Table [T], the total progenitor masses when the collapse 
begins, M^ot, are different from Minit due to the mass loss. The progenitors have the onion- 
like structure: central core is surrounded by shells of lighter elements. We regard the region 
of oxygen depletion as a core hereafter. In Figure |3], the core mass is plotted as a function 
of the initial mass for all models. Since the mass loss rate is larger for high metallicity, the 
models with Z = 0.004 have higher core mass than those with Z = 0.02. The models with 
-^init = 5OM0 become the Wolf-Rayet stars losing large amount of their mass and final core 
masses are somewhat low, which may correspond to type Ic supernovae. As a result, the 
core mass of the model with Minit = SOM© is the highest for each metallicity case. 

To evaluate the spectra of neutrinos emitted from the collapse of above progenitor 
models, which is the goal of this paper, we put the assumptions on the fate of our progenitors. 
For the equation of state by lShen et al. (Il998al jbl). which is adopted below in our simulations. 
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the maximum mass of neutron stars is 2.2Mq in the gravitational mass and 2.5Mq in the 
baryonic mass. We assume that the model with (Minit, Z) = (SOMq, 0.004) forms a black 
hole, because its core mass is larger than the maximum baryon mass. While whether the other 
models can make explosions or not is unclear, we regard them as "supernova" progenitors. 
The correspondence between the initial mass and core mass is not well-established, however, 
the outer region far enough away from the core does not influence the dynamics and neutrino 
emission of the collapse. Therefore, we believe that the systematics discussed in this study 
will be helpful for other progenitor sets if their core masses are known. 



3. Dynamics of core collapse and bounce 



In this section, we follow the early stage of the collapse-driven supernova showing the re- 
sults of our numerical simulation. Here, we utilize the general relativistic neutrino radiation 
hydrodynamics (z/RHD) code, which solves the Boltzmann eq uations for neutrinos together 
with the Lagrangian hydrod ynamics under spherical symmetry (jYamadalll997l : lYamada et al. 
19991 : ISumiyoshi et al.ll2005[ ). We can compute the neutrino transport as well as the dynam- 
ics of collapse. We consider four species of neutrino, z/g, z/g, i^fj. and z/^, assuming that the 
distribution function of z/^ {ur) is equal to that of z/^ (^^)- Neutrino oscillation is not taken 
into account. A detailed description of the numerical simulations such as g eneral relativistic 
hydro dynamics, transport and reaction rates of neutrin os can be found in ISumivoshi et al. 
( 120051 ). Referring the resolution dependences shown in iNakazato et al.l ( 120071 ) . we use 255 
mesh points for the radial Lagrangian coordinate and 20 and 4 mesh points for the energy 
spectrum and angular distribution, respectively. As the initial conditions, the eight progen- 
itor models described in the previous section are adopted an d the outer bounda ry is settled 
far from the core region. We utilize an equation of state bv IShen et al.l (jl998al Jbl). which is 
based on a relativistic mean field model ( jSugahara fc Tokilll994), and inhomogen eous matter 
distribution is described in the Thomas- Fermi approximation (jOyamatsul Il993l ) . 



The density profiles at the time with the central density of 10 g cm" are shown for 
all models in Figure HI We can recognize that their profiles are very similar for the inner 
most IMq and the progenitor dependence is seen in the envelope and outer region of the 
core. The difference will affect the emitted neutrino signal of the collapse- driven supernova. 
As the density increases due to the collapse, the mean free path of neutrinos gets shorter. 
In the collapse stage, the emitted neutrinos are mostly z/g by the electron capture ([2]). Their 
main opacity source is coherent scattering off nuclei; 



Ue + A 



(4) 



When the density exceeds 10 g cm , neutrinos are trapped. The neutrino optical depth 
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d{r) at the radius r is defined as 



d{r) = I 



dr' 



(5) 



J r 



^mfp(r')' 



where Rg is the stellar radius and /mfp is the mean free path of neutrino. Roughly speaking, 
the neutrino trapping occurs inside the neutrino sphere, whose radius, i?^, is defined as 



Once the neutrinos are trapped, the inverse process of the electron capture ([2]) takes place 
and /3-equilibrium for the weak interaction is achieved. As a result, the lepton fraction 
is kept nearly constant and the neutronization is moderated. These trends are shown in 
Figure [5] where the profiles of density, electron fraction and lepton fraction are plotted for 
the model with {Mi^n, Z) = (13Mq,0.02). Here, the lepton fraction Yi is related with the 
electron fraction Ye and the electron-type neutrino fraction Fj,^ as Yi = Ye + Y^^. Note that, 
the electron fraction of the inner region decreases even after the neutrino trapping because 
it is determined by the /3-equilibrium for the given lepton fraction. 

The collapse of the core does not halt until the central density exceeds the nuclear density 
(~ 3 X lO^"^ g cm~^). The infalling core is divided into two parts, to which we refer as the 
inner and outer cores. The inner core contracts subsonically and homologously {v oc r) while 
the outer core infalls supersonically like free-fall (f oc r^^/^). Roughly speaking, it is known 
that the inner core mass corresponds to the Chandrasekhar mass Mch = O.714(yi/O.35)^M0, 
and does not depend on the progenitors. When the inner core density reaches the nuclear 
density, nuclei are closely packed and regarded as uniform matter. Then, the adiabatic 
index increases suddenly due to the repulsive nuclear force and the core restores stability. 
The influence of nuclear repulsion propagates through the inner core and the collapse is 
decelerated. On the other hand, since the pressure wave transmits at the speed of sound, 
the outer core region is still falling supersonically. Therefore, on the boundary of inner and 
outer cores, the shock wave is formed and launched outward. In the meanwhile, the falling 
outer core matter is swept by the shock wave and accretes onto the bounced inner core. This 
compact object is called a proto-neutron star. 

A successful explosion is achieved by the shock breakout through the stellar surface. 
The explosion without the shock stall is called a prompt explosion. However, this scenario 
is not promising according to recent studies and the shock wave is thought to be stalled. 
The shock wave propagation is prevented by the photodisintegration reaction of nucleus ([3]) 
and the neutrino emission. As already mentioned, the photodisintegration reaction is an 
endothermic reaction and consumes the kinetic energy of the shock wave. After the shock 



d{R,) 



2 



(6) 



3 
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wave passes through the neutrino sphere, neutrinos produced by various sorts of reactions 
can escape carrying out the internal energy behind the shock front. Moreover the shock wave 
should overcome the ram pressure of the supersonically infalling outer core. The dynamical 
features described above are seen in Figure M where the snapshots of velocity profile are 
shown for the model with (Minit, Z) = (ISMq, 0.02). 

The stalled shock wave becomes accretion shock and would revive leading to a supernova 
explosion in some way. This is called delayed explosion. Unfortunately, however, the detailed 
scenario is still an open question. Except the QCD mechanism in which a phase transition 
to deconfi ned quark matter na akes the second collapse and shock formation and triggers the 
explosion (IFischer et al.ll201ll ). spherically symmetric (ID) models are thought to fail. In 
this case, the ram pressure of accreting matter is higher and prevents the shock revival. 
On the other hand, if the sphericity is broken, the accretion rate is partially reduced and 
the explosion has an advantage. In fac t, the morphology of collapse-driven supernovae is 
observationally indicated aspherical (e.g. lTanaka et al.ll2009l ). Therefore, most of the recently 
promising scenarios are based on mult i- dimensionality. 

The neutrino heating mechanism has been discussed for years. While, as already men- 
tioned, the neutrinos cool the inner region, they contribute the heating just behind the shock 
front, where free nucleons produced by the photodisintegration reaction ([3]) absorb a small 
part of neutrinos escaping from the inner region. With the aid of this heating, the shock 
wave is considered to revive on a time scale of the order of 100 ms. This scenario was ini- 



tially proposed by iBethe fc WilsonI (119851 ) according to their numerical simulation with ID 
model whereas no other group was able to confirm. Recently, employing multi-dimensional 
simulation, some groups have reported that the onset of the neutrino-driven explosion would 
be helped by aspherical hydrodynamic turbulence such as t he convective instability (e.g. 



Herant et al.lll994j: iFrver fc Warrenll2002l: iMurphv et al. 



2012) and standing accretion shock 



instability (e.g. iBlondin et al.l l2003l : iMarek fc Jankal l2009l : iTakiwaki et all 120121). A nother 
candidate for the explosion scenario is the acoustic mechanism (IBurrows et al.ll2006l ). The 
turbulence generated on the accretion shock travels inward and excites oscillation of the 
proto-neutron star. Then, the acoustic wave originated in the oscillation propagates outward 
and deposits energy on the accretion shock, which causes the shock revival and explosion. In 
this hypothesis, the time scale for the shock revival is estimated as >500 ms, which is longer 
than that of the neut rino heating mechanisrn . Magnetic fields may also be important for 
some supernovae (e.g. iLeBlanc fc Wilsonlll970l : iTakiwaki et al.ll2009t lObergaulinger fc Janka 
20111 ). In fact, neutron stars with strong magnetic fields (magnetars) are observed and their 
progenitors might also have strong magnetic fields. In the MHD mechanism, the rotational 
energy converts to the explosion energy via magnetic fields amplified by the field wrapping 
and/or magnetorotational instability. 
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While, as discussed above, neutrinos would play a key role for a successful explosion, it is 
a small part of them that are absorbed near the shock wave and contribute to the explosion. 
The amount of emitted neutrinos is mainly determined by the released gravitational potential 
of the accreted matter and thermal energy of the nascent proto-neutron star. Since the 
neutrino energy deposition for the explosion has an insignificant effect, our i/RHD simulation 
can be regarded to represent the neutrino signal before the shock revival. In reality, the shock 
wave would move ahead again and the mass accretion stops while when it occurs may depend 
on the explosion mechanism. However, the shock may not be able to revive if the core is too 
massive. It is probable for our model with {Mi^a, Z) = (SOM©, 0.004) and we assume that 
it fails to explode. In this case, we follow the collapse and neutrino emission up to the black 
hole formation. For other models, our z/RHD simulations are terminated at 550 ms after the 
bounce, within which the shock revival is assumed to occur. The evolution after the shock 
revival is dealt in the next section. 

In Figure [TJ we show the evolutions of density, temperature and electron fraction after 
the bounce. They are the results of our z/RHD simulation for the models with (Minit, Z) = 
(13M0,O.O2) and (30Mq, 0.004). For both models, since the proto-neutron star mass gets 
larger by the accretion, the density and temperature rise due to gravitational compres- 
sion. While the density profile is monotonic, the peak of temperature profile resides not 
at the center but at the medium region. This is because the shock wave does not run 
from the center and heats the outer matter, as seen in Figure [61 Since the protons created 
by the photodisintegration reaction cause the electron capture ([2]), the electron frac- 
tion decreases for the shocked region. As recognized in the comparison of the models with 
(Minit, Z) = (I3M0, 0.02) and (BOM©, 0.004), the profiles at the bounce does not depend on 
the progenitor because, as already mentioned, the bounced inner core masses do not differ 
among progenitors. On the other hand, since the accretion rate is determined by the density 
profile of the outer core, the proto-neutron star mass differs among the models for several 
times 100 ms after the bounce. As a result, the proto-neutron star structure is also reflected. 
After the shock propagation, the proto-neutron star settles into a hydrostatic configuration. 
However, the proto-neutron star of the model with (Minit, Z) = (SOM©, 0.004) recoUapses 
suddenly and forms a black hole at 842 n is after the bounce. The qua litative features of 



black hole formation are same as results in ISumiyoshi et al.l ( 120061 . 120071 ) 



4. Proto-neutron star cooling 



We evaluate the neutrino luminosities and spectra in the late phase (after the shock 
revival) by the proto-neutron star cooling (PNSC) simulation except for the black-hole- 
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forming case. In this method, quasi-static evolutions of proto-neutron stars are solved with 
the neutrino transfer b y multigroup flux limited diffusion scheme under spherical symmetry 



with general relativity (ISuzukilll994f ). Here, we follow the hydrostatic structure of the proto- 
neutron star at each time by the Oppenheimer-Volkoff equation while, in z/RHD simulations, 
the equations for hydrodynamics are fully solved. We deal the Boltzmann equations in the 
angle-integrated form for i/g, Ve and v^^ where z/^, z/^, Vr and Vr are treated collectively as v^. 
In contrast, and are treated individually in our z/RHD simulations. Note that and 
have the same type of reactions and the difference in coupling constants is minor. In fact, 
the difference between the distribution functions of and is typically <1% in our z/RHD 
simulations. The equation of state, binning of neutrino energy and neutrino reactions taken 
into account are set to the same with the z/RHD simulation shown in the previous section. In 
this section, we describe the late stage of collapse- driven supernova with our PNSC results. 
Since the time scale of the neutrino diffusion is ~10 s, the evolution is followed till 20 s after 
the bounce. 

We use the results of our z/RHD simulation as initial conditions of our PNSC simulation. 
The central parts up to just ahead of the shock wave are picked up. Here, we take the profiles 
of electron fraction and entropy as functions of the baryon mass coordinate from the z/RHD 
results and, using them, reconstruct hydrostatic configurations with almost steady flow of 
neutrinos, which are used as PNSC initial models. It is confirmed that, except for the close 
vicinity of surface, the obtained density profile is consistent with that of the original z/RHD 
result as expected because the velocity of shocked region is negligible (see Figure^]). Since, as 
already mentioned, detail of the explosion is not known, we set the shock revival time ^revive 
and z/RHD profiles at trevive are used as PNSC initial conditions. Here, we investigate three 
cases as trevive = 100 ms, 200 ms and 300 ms for each progenitor model. When an explosion 
mechanism is assumed, corresponding trevive would be determined. Therefore, we can regard 
that the explosion scenario described above are parameterized by the shock revival time. 
In the PNSC simulation, we follow the evolution of proto-neutron stars without accretion 
because amount of matter which falls back onto the proto-neutron star would be minor after 
the successful explosion. 

In Figure [HI we show the evolutions of density, temperature and electron fraction given 
by our PNSC simulation for the model with (Minit, 2', trevive) = (ISM©, 0.02, 100 ms). One 
can see that the shocked outer mantle of proto-neutron star has relatively high temperature 
and is thermally expanded at initial moment (solid lines). Neutrinos can easily escape from 
there carrying out thermal energy and, therefore, the outer mantle shrinks in a short time 
scale to be denser and hotter while the entropy decreases there (See Figure IH]). Note that 
diffusion fluxes of z/g and v^^ initially transport heat inward from the mantle into the central 
region in parallel outward to the surface because number densities of thermal neutrinos have 
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maxima in the hot mantle. This heat flux attributes to the initial entropy increase in the 
central region. Once the temperature profile becomes monotonic, the entire proto-neutron 
star cools down gradually. In addition to the neutrino cooling, the net flux of electron-type 
neutrinos (z/g — i^e) out of the proto-neutron star carry away the electron-type lepton number 
of the proto-neutron star. This deleptonization just corresponds to the neutronization of 
the proto-neutron star. The nascent proto-neutron star is lepton-rich and composed also 
of considerable protons initially. It becomes the ordinary neutron star with less protons 
achieving the neutrino-less /3 equilibrium in which there remains small amount of protons 
and electrons. The net number flux of electron-type neutrinos is at most ~20% of the number 
flux of Uf.. The decrease of electron fraction shown in Figure [8] represents this process well. 
The density evolution from the onset of collapse to 20 s after the core bounce is illustrated 
in Figure dUl The features of proto-neutron star evolution described above are qualitatively 
common among the models investigated in this study. 

After the revival, the shock wave propagates into the stellar envelope and finally blows it 
off. Since, in the envelope, the energy losses owing to the photodisintegration and neutrino 
emission is quite tiny and the binding energy is subtle, the shock wave is not prevented 
from running outward. While the time scale of the neutrino diffusion is ~10 s, it takes 
from several hours to days for the shock wave to reach the stellar surface. Since supernova 
explosion can be observed optically only after the shock breakout, the neutrino signal reaches 
to the Earth in advance. Finally, the shock wave blasts through the interstellar medium and 
forms a supernova remnant, such as the Crab Nebula. On the other hand, a neutron star 
remains at the center. In Figure [TT], the gravitational mass of neutron stars M^ ns considered 
in our model is plotted. While, according to general relativity, the gravitational mass of a 
neutron star is different from its baryonic mass Mb^NS (see Tabled]), there is a one-to-one 
correspondence between them. We evaluate Mg ns from the ba ryonic mass accreted within 
the shock revival time trevive assuming the equation of state by Shen et al. ( 1998a| jbl). Since 



the accretion rate is higher, ns depends on trevive especially for the progenitors with larger 
initial mass Mjnit. While the neutron star mass could be increased by the accretion from a 
binary companion, the distribution of neutron star masses may give a hint of the explosion 
mechanism. 



5. Neutrino signal 



Neutrinos emitted from the collapse-driven supernova release the gravitational potential 
of the accreted matter and cool the nascent proto-neutron star. For convenience, we divide 
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the neutrino flux after tlie sliock stall, F,y.{E,t), into two terms as 

F,XE,t) = F^\E,t) + F:°"\E,t), (7) 

wfiere subscript i denotes tlie species of neutrinos and E' is a neutrino energy. Tlie first term, 
F^'^{E,t), is an a ccretion term and rela ted with the accretion luminosity, L^^{t), which is 



approximated as ([Thompson et al.ll2003l ) 



GM^{t)M{t) 



where -Rbnd is the outer boundary radius of our z/RHD simulation and Rp{t) is a radius of 
the neutrino sphere at the time t defined by equations ([5]) and (I6l). G, M{t) and Miy(t) 
are the gravitational constant, the mass accretion rate and the mass enclosed by Ru{t), 
respectively. The second term, F^'^°^{E,t), is a cooling term which comes from the thermal 
energy loss of the proto-neutron star. While the accretion term is dominant for the early 
phase (t ~ 100 ms), only the cooling term remains after the shock revival leading the 
explosion. 

As already mentioned, supernova explosion is not successful in most of ID simulations 
owing to the high mass accretion rate. This would be over estimate because some multi- 
dimensional effects such as the convective and standing-accretion-shock instabilities reduce 
the mass accretion in reality. Therefore, we can regard that our ID i/RHD simulation gives 
the maximum case of the mass accretion rate. Thus our z/RHD results of the neutrino flux, 
FI^^^^{E,t), can be regarded as the upper limit: 

= + > F:^^{E,t) + F:^^\E,t) = F,^{E,t), (9) 

where FI^^'^'^^^{E,t) is the maximum possible value of the accretion term. In Figure [12], we 
show the time profiles of luminosities and average energies of emitted neutrinos evaluated 
from i/RHD simulation for all models. The end point of the model with initial mass Mjnit = 
3OM0 and metallicity Z = 0.004 is a moment of the black hole formation. The peak of Ue just 
after the bounce (t = s) corresponds to the neutronization burst. When the shock wave 
propagates through the outer core, nuclei are dissociated into free nucleons and produce a 
large amount of by the electron capture ([2]). A short burst of these neutrinos occurs after 
the shock breakout through the neutrino sphere. This is called neutronization burst. While 
the peak luminosity exceeds 10^^ erg s~^, the duration time is the order of 10 ms and the 
emitted energy is minor comparing with the whole emission of the supernova neutrino. 

The persistent emission after the neutronization burst originates from the mass accre- 
tion and proto-neutron star cooling, as described in equation ([7]). In this stage, u^. and 
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z/g, which are abundantly emitted by the electron and positron captures, respectively, have 
higher luminosity than u^. Nevertheless, the contribution of is not minor because the 
pair processes such as the electron-positron pair annihilation, plasmon decay and nucleon 
bremsstrahlung occur. Thanks to the shock heating, the accreted matter is enough hot for 
electron-positron pair processes. On the other hand, the average energy of is higher than 
those of i^e and Ue- Since does not have charged- current interactions with matter consist- 
ing of no /i^ , , their mean free path is longer at the same position than that of z/g and z/g 
and their neutrino sphere is smaller as recognized by equations and (jH])- Therefore the 
temperature on the neutrino sphere is higher for u^, which makes their average energy also 
higher. As for the progenitor dependence, since the progenitors with higher density for the 
range 1.5-2.OM0 (see also Figure H]) have the higher mass accretion rate, their neutrino lumi- 
nosities are higher as expected from equation (jH]). Their average energies are also somewhat 
higher but the progenitor dependence is not clear especially for the first ~100 ms. Note that 
the neutrino emission stops when the proto-neutron star collapses to a black hole for the 
model with (Mnit,^) = (30Mq, 0.004). 

The results of our PNSC simulation is just corresponding to the cooling term in equation 
([7]). Obviously, they give the lower limit of the neutrino flux because accretion induced 
neutrino flux is not included: 

F™sc(^^i) = F-°'(i?,t) < F^XE,t). (10) 

In Figure [131 we show the time profiles of luminosities and average energies of emitted 
neutrinos evaluated from PNSC simulation for some models. They decrease in time as the 
proto-neutron star cools. The neutrino energy hierarchy {{E^J < {EpJ < (E^,^)) is same as 
that in the accretion phase. The neutrino signals of the models with different initial mass 
Minit and shock revival time trevive are compared in Figure [131 The luminosity and average 
energy are higher for the models with larger neutron star mass M^^ns (see also Figure [TT!) 
while the difference is not so large. 

With inequalities and (ITU]) , we construct the light curve models of neutrino. For 
this, we introduce a ratio factor of the accretion term to its maximum, f{t), as a function 
of time: 

F-'=(E,t) = /(t)F-™(E,t). (11) 

While f{t) may also depend on the species and energy of neutrino, we ignore their depen- 
dences for simplicity. It is required for f{t) to satisfy f{t) ~ 1 for the early phase [t ~ 100 ms) 
and /(t) = for the phase after the explosion. Using f{t), the neutrino flux is expressed as 

F,XE,t) = f{t)F-^-'"^-^{E,t) + F^--\E,t) 

= mF;^^^^iE,t) + il-fit))F^^^^iE,t). 
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The details of explosion dynamics would determine the function f{t). For instance, a neutrino 
signal of the early explosion model corresponds to a rapidly decaying f{t) and a small neutron 
star mass M^ ns- On the other hand, slowly decaying f(t) and large M^ ns give a neutrino 
signal of the late explosion model. 

When supernova neutrinos are actually detected, this study would help us to probe 
the nature of progenitor and remnant. As discussed above, the neutrino luminosity in the 
accretion phase (~100 ms after the neutronization burst) is determined by the progenitor 
model especially for the density profile. The signals in the cooling phase (~10 s after the 
neutronization burst) would provide hints for the mass of remnant neutron star. Moreover, if 
the transition from the accretion phase to the cooling phase is observed, a restriction for the 
explosion mechanism may be possible. Our results can be hopefully utilized as immediately 
comparable templates for a neutrino detection. 

One may be able to use our results for modeling of the supernova neutrino signals. In 
Figure [TU we demonstrate examples of the neutrino light curve and spectrum. They are 
drawn under the assumption, 

1 ; ^ ^ ^revive ~l~ '^shift > 

"S^P I ^, I ) "revive ~r tghift ^ 



'^dccay 

for the model with (Minit, ^revive) = (13M0, 0.02, 100 ms). The junctions of this interpola- 
tion are shown in Figure [151 Since jif) corresponds to the fraction of mass accretion rate to 
its maximum, the decay time scale Tdccay would be a propagation time scale of the revived 
shock wave. When the shock revives at the radius rghock with the escape velocity f esc (''shock), 
it takes 

— 1 /2 

, -Rcore _ p / r shock | -^core A / ^g.NS \ ( shock \ 

tpropagation ~ " ^ ^Gil^;;^ ^ ^ lOOOk^Iy' 

(14) 

for the shock wave to pass the core with the size of -Rcore- Here, Tshock is the most ambiguous 
parameter but at least 100 km < rshock < 1000 km would be satisfied. Therefore, we adopt 
300 km as a typical value and get Xdecay = 30 ms. Note that, this value corresponds to 
also the free-fall time scale of inner core. The time shift tghift is needed for the following two 
reasons. First, Fy^{E, t) is evaluated as the neutrino flux on the outer boundary of our i/RHD 
simulation at the time t. The outer boundary of our PNSC simulation is the proto-neutron 
star surface while that of our i/RHD simulation is located in the envelope. The correction of 
the hght traveling time of the distance between two boundaries (ttravei ~ 0(10) ms) is needed. 
Second, results of PNSC simulations for very early phase are not reliable because, as already 
mentioned, the static density profiles are not consistent with those of i/RHD simulations in 
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the close vicinity of surface. It is needed 0(10) ms for the relaxation. Therefore we set 
^shift = 50 ms, as a typical value. Incidentally, to evaluate the flux numerically, f{t) is settled 
to zero for t > ^revive + tshih + 200 ms, where f{t) ^ 1 in equation f[T51) . Hereafter, for 
simplicity, we take Tdecay = 30 ms and tsmt = 50 ms for all cases. 

In Figure [161 the neutrino number spectra are compared with Fermi-Dirac distributions 
that have the same luminosity and average energy with zero chemical potential. While 
Fermi-Dirac spectra roughly fit our results, the deviation is not small. It is shown that high 
energy neutrinos are also emitted in the accretion phase because the low density outer region 
has high temperature due to shock heating. On the other hand, this region becomes dense 
and the proto-neutron star surface is not heated any more in the cooling phase. Therefore, 
the deficit of high energy neutrinos can be seen in the spectra of late time. Figure [T7] 
compares numerical results and Fermi-Dirac distributions for some cases of time integrated 
neutrino number spectra over the duration of the simulations as in Figure [T6l We find that, 
for the supernova models, the spectrum given by our simulations is relatively well fitted by 
Fermi-Dirac distribution up to ~30 MeV but has high energy tail originated in the accretion 
phase. Even for the black-hole-forming model, the Fermi-Dirac distribution roughly fits our 
spectrum while the numerical uncertainty is larger for the high energy regime. In Figure [TS] 
and Table [H we show the total neutrino energy emitted till 20 s after the bounce for all models 
with the interpolation f|T3l) . Note that, since the model with (Mmit, Z) = (SOM©, 0.004) form 
a black hole, the total neutrino energy emitted till the black hole formation followed by z/RHD 
simulation is plotted. We can see that the total emission energy is related to the core mass 
rather than the initial mass (see also Figure [3]), and it is larger for the explosion models with 
large trevive bccausc the accretion phase is longer and the neutron star mass is larger. 



6. Summary and Discussion 

The purpose of this study is to construct a comprehensive data set of long-term (up 
to ~10 s from the onset of the collapse) supernova neutrino signal for variety of progenitor 
stellar models with different initial masses and metallicities, which would be useful for a wide 
range of research related to supernova neutrinos. To achieve this goal avoiding the difficulty 
of long-term full numerical simulations, we combined two different schemes of numerical 
simulations. The early phase of the collapse-driven supernova, at which the collapsing core 
is bounced and the shock wave is stalled due to the matter accretion, has been followed 
by the general relativistic neutrino radiation hydrodynamics (z/RHD) code. The late phase 
after the shock revival has been dealt by the general relativistic proto-neutron star cooling 
(PNSC) simulation which solves quasi-static evolutions with the neutrino diffusion. The 
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two phases are combined phenomenologically, taking into account the uncertainty about the 
explosion mechanism, and the shock revival time is introduced as a parameter connecting 
the two phases. Although this connection is not perfectly consistent as a single physical 
simulation, this is currently the best way to follow supernova neutrino signals up to ~10 s 
for many progenitor stellar models. There are still many uncertainties about collapse-driven 
supernova physics (e.g., equation of state), but we have chosen the standard or most popular 
parameters, to provide theoretical supernova neutrino emission models expected from the 
standard picture of collapse- driven supernovae. Therefore, the database presented here would 
serve as a standard guideline or template for the supernova neutrino signals. 

It is interestin g to compare our re sult with the past calculations of long-term supernova 



neutrino emission. 



Totani et al.l ( 119981 ) presented the supernova neutrino model till 18 s after 
the bounce for a progenitor with Mjnit = 20Mq and investigated the detectability in detail. 
In their model, the neutrino luminosity decreases by two orders of magnitude within ~10 s, 
which is similar to our results. H owever, their average ener gy gets higher with time, which 



contradicts our results. Recently, iFischer et al.l (120101 12012| ) showed the supernova neutrino 



spectra till 20 s after the bounce for a progenitor with Mi^^t = 18Mq. While they mimicked 
explosion with an artificially enhanced neutrino reaction rate, their results are qualitatively 
consistent with ours. In particular, the drops of neutrino luminosity and average energy due 
to the onset of shock wave revival, where the matter accretion vanishes, are also seen in their 
models. Incidentally, the neutrino average energy is important for the nucleosynthesis such 
as r-process and //-process. 

Because we have calculated many models with different properties of progenitors in a 
consistent manner, we can examine the dependence of supernova neutrino emission properties 
on progenitors. It is quantitatively confirmed that the total emission energy of supernova 
neutrinos is related with the core mass of progenitors. Thus it is larger for the progenitors 
with lower metallicity, but is not monotonically related to the initial mass of progenitors 
due to the mass loss during the pre-coUapse stages. The total neutrino energy emitted also 
depends on the shock revival time that determines the explosion time; it increases with 
^revive because of more material accreting to the collapsed core. The increase is ~ 20-50% 
by changing trevive from 100 to 300 ms. 

In the following we discuss some potential applications of the theoretical supernova neu- 
trino data set for various studies. First, our results could be used as immediately comparable 
templates for the future detection by neutrino detectors. As mentioned repeatedly, the ex- 
plosion mechanism of collapse-driven supernovae is still unknown. The detection of neutrinos 
would give a clue to the diagnosis of the explosion mechanism, because they come from deep 
inside the supernova as shown in this paper. Especially, the drops of luminosity and aver- 
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age energy of neutrinos are important as the observational signature of the explosion. Our 
models would give a quantitative guideline to search such signatures in future detections. 
Furthermore, a sudden stop of neutrino emission would be a signature of a black hole forma- 
tion, like our model with Mi^n = 30Mq and Z = 0.004, providing an exciting opportunity 
to directly observe the birth of a black hole. The neutrino signal reaches to the Earth earlier 
than the electromagnetic signal, and hence triggering the signal by working neutrino detec- 
tors and rapidly informing astronomical communities are of crucial importance. The various 
models presented here would be useful to construct such triggering systems, especially for 
optimizing the efficiency around the detection threshold level. 

It would be possible to predict the spectrum of supernova relic neutrinos using our results 
because we have evaluated the spectra for various progenitor models with different initial 
masses and metallicities. The supernova relic neutrinos are the integration of neutrino flux 
emitted by all collapse-driven supernovae in the causally- reac hable universe. According to 



the s earch for supernova relic neutrinos at Super-Kamiokande (IMalek et al.ll2003l : iBays et al. 



20121). the signal was still not seen but the upper limit was close to the standard predictions 
(IHoriuchi et al.ll2009l ). If the signal of supernova relic neutrinos is actually detected, it would 
give us a unique constraint on the cosmic star formation history and initial mass function. 
A significant excess of diffuse supernova neutrino flux compared with that expected from 
observed supernova rate may indicate a contribution from failed supernovae, most likely 
black hole forming events, like our model with M-mit = SOM© and Z = 0.004 (ILien et al. 
2010h . 



Supernova neutrinos are valuable not only for astrophysics but also for physics of neu- 
trino itself. While we have not taken into account the neutrino oscillation in calculations of 
the models presented here, it can easily be de alt as a post-process to predict final neutrino 
signal reaching detectors on the Earth, (e.g. iKotake et al.l 120061 ). Recently, as indicate d 
by the results from T2K (lAbe et al.l 120111 ) and MINOS experiments (lAdamson et al.l |20 111 ), 
the mixing angle of neutrino oscillation 613 is confirmed to be nonzero and eva luated as 
sin^6'i3 ~ 0.1 by the results fro m reactor neutrino e xperiments such as Daya Bay (lAn et al. 



2OI2I), RENO flAhn et al.ll2012l ) and Double Chooz flAbe et al.ll2012l ). Thus, at present, the 
most undetermined parameter in the neutrino oscillation is the mass hierarchy. Since, for 
sin^26'i3 ~ 0.1, the survival probabilities of Ue, and Vr in the stellar envelope are different 
for the normal and inverted mass hierarchies, det ections of the super nova neutrinos would 
give useful information of the mass hierarchy (e.g. iKotake et al.l 120061 ). 



The supernov a neutrino signal depen ds on the equation of state, especially for the black- 



hole-forming case ( ISumiyoshi et al. 



20061). The effect of hyperons, which is not taken into 



account in our equation of state (IShen et al.l ll998al Jbl). would be important not only for 
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the black hole forma tion ( iNakazato et al.ll2012l ) but also for the proto-neutron star cooling 
( iKeil fc Jankal Il995l ). However, hyperonic equation of state is an unsettled hot topic. In 
particular, recentl y, the mass of the bin ary millisecond pulsar J1614-2230 was evaluated 
as 1.97 ± O.O4M0 (jPemorest et al.l 120101 ). Unfortunately, this remarkable precision thanks 
to a strong Shapiro delay signature excludes almost all models of hyperonic equation of 
state, because the rnaximum mass of n eutron stars gets lower by the hyperon inclusion (e.g. 



Ishizuka et al.l l2008t IShen et al.l 120111 ) . Neutrino interactions in matter are also affected 



by equation of sta te. In this study, we have adopted a single nuclear equation of state of 



Shen et al.l (ll998al Jbl) without hyperons, which is based on mostly standard assumptions and 



can be regarded as a baseline model. While there are some issues beyond the scope of this 
study, we hope that the result of this paper will be useful for further progress of the related 
fields in astrophysics and neutrino physics. 
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A. Supernova Neutrino Database 

The numerical data of supernova neutrino emission computed in this study are publicly 
available on the Web at 



|http : / / asphwww . ph . noda . tus .ac.jp/ snn/ 



This data set is open for general use in any research for astronomy, astrophysics, and physics. 
Not only the original data of i/RHD and PNSC simulations but also combined data from the 
onset of collapse to 20 s after the core bounce with the interpolation (IT^ are provided. On 

the Web, the differential neutrino number flux ^''^^^ and differential neutrino number 

luminosity ^''^'^ ""^ at the time tn are prepared for each neutrino species Ui and energy 

AEk 
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bin Ek- Moreover, the spectral data integrated from the onset of collapse to 20 s after the 
core bounce with the interpolation f|T3|) are also shown. 
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Note. — Minit and Z are the initial mass and metallicity of progenitors, respectively. Mtot, Mhb and Mco are a total progenitor mass. He core mass and CO 
core mass when the collapse begins, respectively. Since models with Minit = 5OM0 become Wolf-Rayet stars, Mne is not defined and Mco equals Mtot. Moore is a 
core mass which is defined as the region of oxygen depletion, trevive is a shock revival time. M^^ns and Mg ns are a baryonic mass and gravitational mass of the 
remnant neutron stats, respectively. The mean energy of emitted Vi till 20 s after the bounce is denoted as (^Ei,^) = Eiy.^tot/N„.^tot, where -Ei/.,tot and N„.^tot are 
the total energy and number of neutrinos. Ux stands for /i- and r-ncutrinos and their anti-particles: E^^ = E^^ = Et?^^ = E,^^ = Ei^^. -Ei/.^^.tot is a total neutrino 
energy summed over all species. The model with Minit = 3OM0 and Z = 0.004 is a black-hole-forming model, for which mean and total neutrino energies emitted 
till the bla<;k hole formation are shown. 
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Fig. 1. — Density profiles of progenitor models with the metallicity Z = 0.02 (left panel) 
and Z = 0.004 (right panel). In both panels, upper left, upper right, lower left and lower 
right plots corresponds to the models with the initial mass Minit = 13Mq, 20Mq, 30Mq and 
50Mq, respectively. Dotted vertical lines are shown to guide eyes. 
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Fig. 2. — Composition profiles of progenitor models with the initial mass Minit = ISM© and 
the metallicity Z = 0.02 (left panel) and Mnit = 30Mq and Z = 0.004 (right panel). Plots 
are shown as a function of the baryon mass coordinate (enclosed mass). In both panels, 
">Si" means a sum for the elements heavier than silicon. 



- 26 - 



2.8 







2.6 - 



2.4 



m 2.2 

m 

03 



CD 

o 

o 



2 - 



1.8 - 



1.6 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




Z = 0.02 



I ^ I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

10 20 30 40 50 60 
initial mass (Mg) 



Fig. 3. — Core mass of progenitor models as a function of the initial mass. Solid and 
dashed lines corresponds to the cases with the metallicity Z = 0.02 and 0.004, respectively. 
The horizontal dott ed line represents t he maximum baryonic mass of neutron stars for the 
equation of state by IShen et al. 
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Fig. 4. — Density profiles at the times with the central density of 10^^ g cm~^ for progenitor 
models with the metallicity Z = 0.02 (left panel) and 0.004 (right panel). In both panels, 
sohd, dashed, dotted and dot-dashed hnes correspond to the models with the initial mass 
Minit = 13M0, 2OM0, 3OM0 and SOM©, respectively. 
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Fig. 5. — Snapshots of the collapsing core for the model with the initial mass Minit = ISM© 
and the metallicity Z = 0.02. The left, center and right panels show the density, electron 
fraction and lepton fraction profiles, respectively. In all panels, solid, dashed, dotted and 
dot-dashed lines correspond to the times with the central density of 10^^ g cm~^, 10^^ g cm~^, 
]^q13 g f.T^~3 ^T^^ ]^gi4 g Q]2i-3^ respectively. 
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Fig. 6. — Snapshots of velocity profile for the model with the initial mass Mjnit = ISMq and 
the metallicity Z = 0.02. The labels represent the time ordering, where 4 and 8 correspond 
to the time at the bounce and 5 ms after the bounce, respectively. 
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Fig. 7. — Evolutions of the density, temperature and electron fraction profiles by the simu- 
lation of neutrino radiation hydrodynamics. The upper plots show the results of the model 
with the initial mass Minit = 13Mq and the metallicity Z = 0.02, and the lower plots do 
those of the model with Minit = SOMq and Z = 0.004. In all panels, solid, dashed and dotted 
lines correspond to the times at the bounce, 100 ms after the bounce and 300 ms after the 
bounce, respectively. The dot-dashed lines represent the profiles at 550 ms after the bounce 
and the profiles at 840 ms after the bounce (2 ms before the black hole formation) for upper 
and lower panels, respectively. 
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Fig. 8. — Evolutions of the density, temperature and electron fraction profiles by the simu- 
lation of proto- neutron star cooling for the model with the initial mass Mjnit = 13Mq, the 
metallicity Z = 0.02 and the shock revival time ^revive = 100 ms. In all panels, sohd, dashed, 
dotted and dot-dashed lines correspond to the times at 100 ms, 1 s, 7 s and 20 s after the 
bounce, respectively. 
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Fig. 9. — Snapshots of entropy profile for the model with the initial mass Minit = ISMq, the 
metallicity Z = 0.02 and the shock revival time trevive = 100 ms. The notations of lines is 
the same as Figure |8l 
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Fig. 10. — Evolution of the density profile for the model with the initial mass Mjnit = 13Mq, 
the metallicity Z = 0.02 and the shock revival time ^revive = 100 ms. Solid, dot-dashed and 
dashed line correspond to the onset of collapse, 100 ms after the bounce and 20 s after the 
bounce, respectively. Dotted vertical lines and horizontal arrows are shown to guide eyes. 
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Fig. 11. — Gravitational mass of neutron star models with the metallicity Z = 0.02 (left 
panel) and 0.004 (right panel) as a function of the shock revival time. The notations of lines 
is the same as Figure H] in both panels. 
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Fig. 12. — Luminosities (upper plots) and average energies (lower plots) of the emitted neu- 
trinos as a function of time after bounce from the z/RHD simulations. The panels correspond, 
from left to right, to z/g, z/g and (= ^^j.-, ^t)- The results for the models with the 

metallicity Z = 0.02 are shown in the top panels, and those for the models with Z = 0.004 
are shown in the bottom panels. In all panels, solid, dashed, dotted and dot-dashed lines 
correspond to the models with the initial mass Mjnit = 13Mq, 20Mq, SOMq and 50Mq, 
respectively. "BH" means a black-hole-forming model with Mjnit = 30Mq and Z = 0.004 
and its end point corresponds to the moment of black hole formation. 
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Fig. 13. — Same as Figure [12] but from the PNSC simulations. In the left panel, signals of 
Vf, (solid lines), Vf. (dashed lines) and z/^. (dot-dashed lines) are shown for the model with 
(Minit, Z, ^revive) = (13Mq, 0.02, 100 ms). In the central panel, v^. signal are shown for the 
models with (Z, trcvive) = (0.02, 100 ms) and Minit = 13Mq (solid lines), 2OM0 (dashed lines), 
3OM0 (dotted lines) and SOM© (dot-dashed lines). In the right panel, z/g signal are shown for 
the models with (Mjnit, 2') = (I3M0, 0.02) and trevive = 100 ms (sohd lines), 200 ms (dashed 
lines) and 300 ms (dot-dashed lines). 
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Fig. 14. — Time evolution of neutrino luminosity and average energy (left), and number 
spectrum of (right) from i/RHD and PNSC simulations with the interpolation (fT3|) for 
the model with (Mjnit, Z, trevive) = (ISMq, 0.02, 100 ms). In the left panel, solid, dashed and 
dot-dashed lines represent u^, and (dot-dashed lines), respectively. In the right panel, 
the lines correspond, from top to bottom, to 0.1, 0.25, 0.5, 2, 4 and 15 s after the bounce. 
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Fig. 15. — Time evolutions of neutrino luminosity and average energy of z/g for the model 
with (Minit, -Z') = (13M0,O.O2). Thick dashed and thick dot-dashed lines represent the 
interpolations ( fT3|) with trevive = 100 ms and trevive = 300 ms, respectively. Thin solid lines 
show the results of z/RHD and thin dashed and thin dot-dashed lines do the results of PNSC 
simulations with trevive = 100 ms and trevivc = 300 ms, respectively. 
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Fig. 16. — Number spectra of at selected times for the model with (Mjnit, 2', trevive) = 
(13M0, 0.02, 100 ms). Solid, dashed and dot-dashed lines correspond to 0.1, 0.5 and 15 s 
after the bounce, respectively. Thick lines show the results of our simulations while thin 
lines are Fermi-Dirac spectra with the same luminosity and average energy as the numerical 
results. The chemical potential is set to zero for the Fermi-Dirac distribution. 
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Fig. 17. — Time integrated number spectra of Ve over the duration of the simulations. The 
left, center and right panels correspond to the models of supernova with (Mmit, trevive) = 
(13Mq, 0.02, 100 ms), supernova with (Mi^it, Z, t^evive) = (3OM0, 0.02, 300 ms) and black hole 
formation with (Minit, Z) = (3OM0, 0.004), respectively. The meaning of thick and thin lines 
are the same as Figure fT6l 
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Fig. 18. — Total neutrino energy emitted till 20 s after the bounce for the models with 
the shock revival time trevive = 100 ms (left), 200 ms (center) and 300 ms (right). They 
are computed from z/RHD and PNSC simulations with the interpolation ( fT3ll except for the 
model with the initial mass Minit = 30Mq and the metallicity Z = 0.004, for which the 
neutrino emission till the black hole formation followed by z/RHD simulation is plotted in all 
panels. The notations of lines is the same as Figure [31 



